Optimal linear optical implementation of a single-qubit damping channel 
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We experimentally demonstrate a single-qubit decohering quantum channel using linear optics. 
We implement the channel, whose special cases include both the amplitude-damping channel and the 
bit-flip channel, using a single, static optical setup. Following a recent theoretical result [M. Piani 
et al, Phys. Rev. A, 84, 032304 (2011)], we realize the channel in an optimal way, maximizing 
the probability of success, i.e., the probability for the photonic qubit to remain in its encoding. 
Using a two-photon entangled resource, we characterize the channel using ancilla-assisted process 
tomography and find average process fidelities of 0.9808 =b 0.0002 and 0.9762 ± 0.0002 for amplitude- 
damping and the bit-flip case, respectively. 

PACS numbers: 42.50.-p, 42.50.Ex, 03.67.-a, 03.65.Yz 



Introduction. Time evolution in quantum mechanics 
converts a density matrix to another density matrix. This 
evolution is referred to as a quantum channel and can be 
described mathematically as a completely positive (CP) 
map [1] . Because of the generality of the concept of quan- 
tum channels, their use is ubiquitous in quantum infor- 
mation. For example, unitary quantum channels are used 
in quantum computing to describe quantum gates. Non- 
unitary channels, on the other hand, describe the inter- 
action of quantum states with an environment, and have 
recently been connected to fundamental physical ques- 
tions in quantum information science, such as channel 
capacity, superadditivity [2j [3] and bound entanglement 

Linear optics and single photons have several charac- 
teristics that make them an ideal testbed for quantum 
information. Single-qubit unitaries are easy to imple- 
ment as, for polarization encoded qubits, they only re- 
quire waveplates. Photonic qubits also exhibit long co- 
herence times, and spontaneous parametric down con- 
version allows the generation of high-quality entangled 
states, which can be easily manipulated. However, cer- 
tain operations, such as the two-qubit CNOT-gate, are 
difficult in this architecture |5| |6^ , and can only be im- 
plemented probabilistically [7H9l. 

Unfortunately, the ease of single-qubit operations does 
not extend to more general CP maps. Some quantum 
channels, like the depolarizing single-qubit channel [1] 
can be implemented with unit probability, but this is not 
the case in general. For instance, the amplitude-damping 
channel, a non-unital quantum process, has been imple- 
mented in linear optics only with a limited success proba- 
bility of 1/2 jTOl[TT]. In one experiment [12 , the interac- 
tion of a qubit with ancilla modes was implemented such 
that the resulting counts replicated those expected for 
an amplitude-damping channel. However, the quantum 
information did not remain in its original encoding and 
cannot be viewed as a single-qubit channel. Recently, 
it was shown by Piani et al [13 that, any single-qubit 



quantum channel could be implemented probabilistically 
using linear optics and postselection, i.e., similar to many 
two-qubit operations. Moreover, they proved that such 
implementations can be achieved with the optimal suc- 
cess probability. 

In the present work, we use this recent theoretical re- 
sult to design and demonstrate a linear-optics-based im- 
plementation of a certain class of non-unital single-qubit 
quantum channels called "damping channels" . The class 
of channels we focus on can be parametrized by two real 
numbers: a, f3. In the operator-sum representation, the 
channel's action on an arbitrary quantum state p can be 
written as E{p) = AipA\^ where the two Kraus oper- 
ators are p!4] : 

. ( cos a \ . / sin/3\ 
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This channel is of great interest as its special cases in- 
clude the amplitude-damping (a = 0) and bit-flip (a = 
/3) channels, both of which are common sources of error 
in other implementations of quantum information pro- 
cessing, such as ion traps. Furthermore, it belongs to 
the small class of quantum channels for which the quan- 
tum capacity can be directly calculated via the coherent 
information [15 . 

Here, we experimentally realize this single-qubit damp- 
ing quantum channel using linear optics. We can use the 
setup to add controlled amounts of noise of various types 
to a single qubit. The schematics of the experimental 
setup are shown in Fig. [T] The key step in the imple- 
mentation of the channel is the splitting of the polar- 
ization encoded information into different spatial modes, 
which then allows for the manipulation of different logical 
states independently of one another. An arrangement of 
half-wave plates and liquid-crystal retarders allows us to 
probabilistically implement both Kraus operators within 
a single, static optical setup. We characterize the chan- 
nel using a new ancilla-assisted quantum process tomog- 
raphy method, and show the optimality of our optical 
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FIG. 1: The experimental setup. We use spontaneous 
parametric down conversion in periodically poled KTiOP04 
(PPKTP) to generate entangled photon pairs of the form 
1$^) = (\HH) + I W^)) which are subsequently coupled 
into single-mode fibres. One of the photons is sent through 
the damping channel parameterized by a and /3, which are set 
by the angles a, 6, c and d of four half- wave plates (HWPs). 
Two liquid-crystal retarders (LCRs) switch ant i- cor relatively 
between the identity, 1, and the Pauli X operation. The po- 
larization of each photon is measured by an analyzer (A and 
B) consisting of a half- and a quarter-wave plate (QWPs) fol- 
lowed by a polarizing beam splitter (PBS). Eventually, the 
photons are detected by single-photon counting modules. 



implementation, with our success rates in the amplitude- 
damping case surpassing those of previous implementa- 
tions [10 , 11 . In order to characterize the action of the 
channel on entanglement we study the amount of entan- 
glement of photon pairs when one photon is sent through 
the channel. 

Optimality of the implementation. Following Ref. [13], 
it can be shown that the probability of success for 
a specific Kraus decomposition {Ai} is Psucc({^i}) = 

(Zlill^illL) ' where the norm ||M||^ is the largest 
singular value of the operator M. Maximizing over all 
possible Kraus decompositions Ai describing the chan- 
nel allows to achieve the optimal success probability 
Psucc = max^. ^ ||V||2 • For our particular channel, if 

II * II c>o 

we assume that cos(a) > cos(/3), this yields: 
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In order to achieve Pgucc h.dive to implement each 
Kraus operator with individual probabilities pAi = 
ll^dlL ■ -Psucc- We find that the optimal probability of 
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success is achieved for pAn = — cos a _ 

■ 2 o 

cos^^o+sin^ 13 ' show experimentally that for various 

values of a and /3, which can be independently controlled 
in our experiment, we indeed achieve this upper bound. 

Ancilla- assisted quantum process tomography. Quan- 
tum process tomography (QPT) allows to experimentally 
reconstruct the superoperator describing an unknown 
physical process. Ancilla-assisted QPT (AAQPT) uses 



ancillary qubits to facilitate the reconstruction procedure 
for quantum state measurements. It has also been shown 
[16] that AAQPT gives decreasing statistical errors as the 
amount of entanglement between the primary and ancilla 
systems is increased. 

AAQPT has been used to study various unitary quan- 
tum gates [16] but has not yet been extended to the char- 
acterization of non-unital channels. In our work, and in 
contrast to previous AAQPT schemes [TTHTQ] , we do not 
assume an ideal description of our initial state, but rather 
measure and include it when using a maximum-likelihood 
technique to find the physical matrix that best describes 
the action of the experimentally implemented channel. 

The standard techniques for QPT and AAQPT are 
described in [1 and [16], respectively. Below, we out- 
line our method following their nomenclature. Consider 
a two-qubit state, Pab^ whose density matrix is known; 
e.g., it might have been reconstructed using quantum 
state tomography (QST). The quantum channel £ acts 
on subsystem A, while subsystem B is unaffected. The 
transformed two-qubit state after the channel is = 

(g) 1) (pab)- Characterizing Pab^ ^-S-? by performing 
standard QST, allows for reconstruction of the process 
using the Choi-Jamiolkowski isomorphism [20l [21] . 

The quantum process can then be written as = 

El,nLo Xmn{Em ^ l)pAs(^n ^ whcrc {Ei} are op- 
erators which form a basis in the space of d x d matrices 
= 2 in our case). It is common to use the basis formed 
by the Pauli matrices y, Z}. The -dimensional 

process matrix x then fully describes the quantum pro- 
cess. In our maximum-likelihood technique, we parame- 
terize X by — (i^ = 16 real numbers [181 IE] and seek 
to minimize the following function: 



i=l 



(3) 



such that the resulting x most closely resembles a phys- 
ical quantum process. Here, i labels the measurement 
setting in the final QST, u is the number measurement 
settings, Ui is the number of two- fold coincidence counts 
recorded in the i^^ setting. A/" corresponds to the number 
of photons incident on the detectors. Mi is the projector 
of the i^^ measurement, and A is a Lagrange multiplier 
used to force the resulting process matrix to be trace 
preserving [17 . 

Experiment. We use the experimental setup shown in 
Fig. [T] to implement the quantum channel defined by the 
Kraus operators in Eq. [ij Two 40 mm calcite beam dis- 
placers are used to construct an interferometer. Within 
these beam displacers, photons with horizontal {\H)) and 
vertical {\V)) polarization are spatially displaced with re- 
spect to each other [9 . Half- wave plates (HWPs) are 
used to set the amount of damping by allowing to adjust 
a and P in Eq. [ij The relations between these parame- 
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FIG. 2: (a) Probability of success as a function of damping 
parameter /3 for cases a = 0, a = |/3 and a — ^. The shaded 
region below 0.5 represents probabilities of success in previous 
optical implementations of the amplitude-damping channel, 
see Refs. [10 -12 (b) The tangle, r, of the resulting two-photon 
state as a function of damping ^ after one photon has passed 
through the damping channel. Errors in the experimental 
data are calculated from Poissonian noise in the coincidence 
counts and are not visible on the scale of the plots. The 
solid lines in both panels represent the theoretically expected 
dependance. In (b), the experimental density matrix of the 
input state was used for the calculations. The shaded regions 
around the theory curves in both plots represent the expected 
standard deviation assuming 1° and 1% rotation errors in the 
HWPs and LCRs respectively. The margin of error in the 
a — P case is significantly smaller than for the other cases 
due to the fact that the HWP angles and LCR settings all 
lie at points where partial derivatives of the Kraus-operator 
matrix elements are zero. 



ters and the individual angles a, 6, c, d of the four HWPs 
are given by sin 4a = 5 = ^ _ ^ sin 4c = — 

" ^ cos ol' 4 " sm p 

and d = f — c. The channel is realized by switching 
randomly between the Kraus operators Aq and Ai. The 
switching is performed using two liquid-crystal retarders 
(LCRs). We set the two LCRs to Xi and 12, respec- 
tively, to implement Aq^ and we set them to li and X2 
to implement Ai. Here, the subscripts represent the ac- 
tion of the first and the second LCR. The probabilities, 
PAq and p^^, with which each configuration is realized 
are determined by the values of a and {3 such that the 
overall success probability of realizing the channel is op- 



timal [13]. The switching rate of the LCRs was chosen 
to be 10 Hz, significantly faster than the integration time 
for a single measurement (5 s). 

To characterize the channel, we use the AAQPT 
scheme outlined above. Our resource state is an entan- 
gled photon pair generated in a type-II down conversion 
source in a Sagnac configuration [22l|23]. A 0.5 mW laser 
diode at 404.5 nm pumps a 25 mm periodically-poled 
crystal of KTiOP04 (PPKTP). This typically yielded a 
coincidence rate of 10 kHz. The characterization of the 
channel is executed as follows: The HWP angles a, 6, c 
and d are set to zero and the LCRs to Xi and I2 such 
that the channel acts as the identity map. A Q ST is per- 
formed to obtain the density matrix of the input state, 
pAB' The HWP angles and the probabilities for switch- 
ing the LCRs and the HWP angles are then set accord- 
ing to the values of a and /3. Another QST yields the 
output state, p'ab- QST involves recording coincidences 
for all combinations of the eigenstates of the Pauli X, Y 
and Z operators. For each of these 36 projective mea- 
surements, we integrated coincidence counts for 5 s. The 
resulting data were then used in conjunction with Eq. [3] 
to reconstruct the superoperator describing the quantum 
process. 

Results. We now turn to our main result, the opti- 
mality of our quantum channel implementation. Fig. [2^ 
shows the probability of success for the amplitude- 
damping, bit-flip, and one in-between case {a = |/3). 
Since amplitude-damping manifests itself as photon loss 
in our particular implementation, determining the prob- 
ability of success reduces to measuring the transmission 
of the channel. The experimental data closely follow 
the theoretical predictions (solid lines) that are based 
on Eq. (2] 

Previous optical implementations of the amplitude- 
damping channel [lOl [TT] have given at most 50% prob- 
ability of success [24], whereas here we find that only in 
the case of maximum damping (/3 = 7r/2) the probability 
of success decreases to 50%. The experimental results 
for the success probability closely resemble the theoret- 
ical prediction. This is also true for our experimental 
implementation of the bit-flip channel and the a = |/3 
case of single-qubit damping. 

Fig. shows the tangle [25 of the two-photon out- 
put density matrix, p^^, for the amplitude-damping, 
bit-flip, and a = |/3 cases, where one of the two pho- 
tons passes through the quantum channel. Theoretical 
curves are based on the action of the respective ideal 
quantum channels on the experimental input state. It 
can be seen in each case, that the experimental data 
agrees well with the theoretical prediction, showing that 
the experimentally implemented channel closely resem- 
bles the ideal one. The process fidelity is defined by = 
T^'^yQT^^Xid^/T^ [IS], where Xexp and Xid are the ex- 
perimental and ideal process matrices, respectively. We 
find process fidelities of 0.9808 ± 0.0002, 0.9762 ± 0.0002, 
0.9805 ± 0.0003 for the a = 0, a = /3 and a = |/3 cases 
of the channel, respectively. We also compute the so- 



4 




Re(Xexp) Re(Xid) 



quantum channel with hnear optics in which both type 
and amount of decohering noise can be precisely con- 
trolled. A single, static optical setup can perform as the 
amplitude-damping channel, the bit-flip channel, or more 
general cases characterized by two real parameters, a and 
p. Most importantly, we have shown that this channel 
has been implemented in an optimal way, so as to maxi- 
mize the probability of success. The channels were char- 
acterized using a new approach to ancilla-assisted process 
tomography and, in all cases, operate with high fidelity. 
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FIG. 3: (a) Measured process fidelity (solid data points) and 
trace distance (unfilled data points) as a function of damping 
for each of the three cases studied. Error bars (^ lO""^), 
calculated using Monte-Carlo simulations adding Poissonian 
noise to the measured state tomography counts in each run, 
are too small to see on this scale, (b) Real and imaginary parts 
of the experimentally determined and ideal process matrices 
at maximum amplitude-damping (a = 0, /3 = 7r/2). 

called maximum trace distance [T, which is defined as 
V = maxp.^ ^Tr Ip^^p - pjf^tl^ where \A\ = VA'^A and 

Pout^'"^ = Em,n Xmn^'^^^m Pin^i- Operationally, V cor- 
responds to the highest probability of distinguishing be- 
tween the experimental and ideal channels using the best 
possible input state. The average maximum trace dis- 
tance over all damping values measured were found to 
be 0.1028 ± 0.0009, 0.1077 ± 0.0006, 0.1075 ± 0.0009 for 
the three respective cases. In Fig. |3] we show both the 
process fidelity and maximum trace distance as a func- 
tion of damping for each case, as well as the reconstructed 
X matrix for a = 0, /3 = 7r/2. 

Summary. Decoherence plays an important role in 
quantum information science. Investigating its effects 
requires careful and well-controlled implementations of 
these noisy processes. Non-unital damping channels, like 
the ones studied here, are crucial in further understand- 
ing quantum communication, in determining channel ca- 
pacities and for the generation of bound-entangled states. 
We have implemented a general damping single-qubit 
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